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Abstract 

We discuss small perturbations on the self-accelerated solution of the DGP 
model, and argue that claims of instability of the solution that are based on 
linearized calculations are unwarranted because of the following: (1) Small 
perturbations of an empty self-accelerated background can be quantized con- 
sistently without yielding ghosts. (2) Conformal sources, such as radiation, do 
not give rise to instabilities either. (3) A typical non-conformal source could 
introduce ghosts in the linearized approximation and become unstable, how- 
ever, it also invalidates the approximation itself. Such a source creates a halo 
of variable curvature that locally dominates over the self-accelerated back- 
ground and extends over a domain in which the linearization breaks down. 
Perturbations that are valid outside the halo may not continue inside, as it 
is suggested by some non-perturbative solutions. (4) In the Euclidean con- 
tinuation of the theory, with arbitrary sources, we derive certain constraints 
imposed by the second order equations on first order perturbations, thus re- 
stricting the linearized solutions that could be continued into the full nonlinear 
theory. Naive linearized solutions fail to satisfy the above constraints. (5) Fi- 
nally, we clarify in detail subtleties associated with the boundary conditions 
and analytic properties of the Green's functions. 
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1 Introduction and summary 



The accelerated expansion of the Universe could be due to modification of gravity 
at cosmological distance scales. One "existence proof" for such a scenario is the 
DGP model in which the self-accelerated solution can be used to describe the 
observed speed-up of the cosmic expansion P]. 

While there are a number of interesting works exploring the consistency of the 
self-accelerated solution with the data, see, e.g. [3]-[Tn], the linearized perturbations 
on this background were claimed to contain a ghost- like mode fK^ 1^ 1^ 1^ 1^. 

If the results of the linearized calculations were reliable, this would exclude the 
self-accelerated background from a set of physically acceptable solutions of the DGP 
model. However, we will show below that the issue is more involved. Small pertur- 
bations on an empty self-accelerated background can be quantized without yielding 
ghosts. Moreover, for conformal sources, such as radiation, the small perturbations 
are still ghost free. Non-conformal sources, on the other hand, could introduce ghost- 
like instabilities in the linearized theory. However, in the latter case the results of 
the linearized approximation are unwarranted; no conclusion on the (in) stability of 
the solution of the DGP model can be drawn without invoking non-linear dynamics. 

An important role of non-linearities in this model have been emphasized in a 
somewhat different context. In a simplest case of gravity of a static source, e.g., 
the Sun, it is the breakdown of naive linearization that enables one to evade [23] 
the van Dam-Veltman-Zakharov (vDVZ) discontinuity which, if present, would 
have ruled out the model already by the Solar system data, without any further 
need to consider cosmology. 

Fortunately, the breakdown of the linearized perturbation theory does not neces- 
sarily mean incalculability. Instead, an expansion in other small parameters should 
be carried out, see, e.g., j2Sl , EZ] , EH] • For instance, despite the fact that the nearby 
gravitational field produced by a static source (such as Sun, Earth, etc.) is incalcu- 
lable within the conventional linearized approximation, the alternative approaches 
give rise to a weak field which approximates the result of General Relativity (GR) 
to a high degree of accuracy, with tiny deviations that could be potentially mea- 
surable (see E] and ^^l)- The fundamental reason for this is the absence of 
a continuous transition from DGP to GR in the linearized theory (i.e., the vDVZ 
discontinuity is present there), whereas, in a full non- linear model such a transition 
does exist, at least in exact solutions of the classical theory 

After these general remarks, we turn below to a brief description of the content 
and main results of the present work, with all the important technical details left 
for the bulk of the paper and Appendices. 

We start in Section 2 by studying an empty self-accelerated solution, i.e., the 
background on which no external matter/radiation or non-linear gravitational sources 
are introduced. This is an artificial reduction of the theory, but it serves as an 
instructive playground. We look at small gravitational perturbations on this back- 
ground. The linearized approximation gives a well-defined spectrum of tensor per- 
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turbations, but in addition, there is a scalar mode that could or could not give 
rise to the ghost-like instability. We will show that for the empty brane the scalar 
acts as a Lagrange multiplier. Such a theory is analogous to conventional QED in 
which the Lorentz gauge-fixing condition is introduced by the Lagrange multiplier 
technique. It has been known that such a model can be quantized consistently, 
e.g., in the canonical approach by the Nakanishi-Lautrup method [S^EDIj ot more 
straightforwardly, by the path integral approach. Likewise, in the gravitational case 
at hand, we will show that the perturbations can be quantized in a path integral 
approach without negative norm states (i.e., without the ghost). We will discuss 
why our results differ from those conclusions of Ref. j2S| that deal with the empty 
self-accelerated solution. 

Things change when non-conformal matter or non-linear gravitational sources 
are allowed in the theory. For simplicity we consider below an external matter 
source Tfj^u- The linearized calculations suggest that T^j, with a nonzero trace could 
excite a ghost-like instability^. If this were all, then any nonzero T could start to 
radiate the ghost-like mode, rapidly destabilizing itself. 

The above logic, however, misses the main feature of the DGP model, the im- 
portance of the nonlinear interactions. To put our discussions in a context let us 
first consider an empty brane with fiat space-time as a background (this is not a 
self-accelerated solution yet!). Moreover, consider for simplicity a spherically sym- 
metric matter distribution of radius tq and density po- The 4D gravitational radius 
of this source is vm = 2GnM = 2G'Ar(4/3)7rroPo- There is a new distance scale in 
the theory [22] (called the Vainshtein scale 

U^{rM/Hiy'\ (1) 

where Hq ~ 10^^^ GeV is the present-day value of the Hubble parameter. Below 
the scale (P), i.e., for r ^ r* the linearized approximation breaks down For 
this breakdown to be relevant, has to be bigger that tq. This translates into 
the condition on the density of the source po ^ where pc ~ IQ'"^^ g / cm? is a 
critical density of the Universe. For such sources the invariant scalar curvature, 
i?, is bigger than the space-time curvature of the observable Universe, i.e., R > 
Hi. Most of the sources around us satisfy this condition; we will refer to them as 
supercritical sovLTces. In any realistic setup, and on the self-accelerated background in 
particular, the scalar curvature due to the background is already "critical" R ~ H^; 
therefore, within any overdensity on top of this background, the local curvature 
will be supercritical R > H^, and one should expect the naive linearization to break 
down. One could certainly think of a case when one "prepares by hand" a subcritical 
source by distorting locally the self-accelerated background. Such a source, however, 
would by itself be unstable. One way for this instability to develop, is for the 

■^This statement, even within the Unearized approach itself, needs further clarification; for details 
see calculations and discussions in Section 4. 

similar scale was first found by Vainshtein in the context of 4D massive gravity |31 | 
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source to reduce its size until its density becomes supercritical. This is energetically 
favorable since the bare mass of the supercritical sources gets screened |13j: A 
supercritical source responds to the breakdown of the linearized approximation by 
creating a halo of a variable curvature which extends over a domain or size r*. This 
halo screens the stress-tensor of the source. Moreover, the local curvature in the 
halo dominates over the curvature of the self-accelerated background itself |13j . 

How can nonlinear-effects change the conclusions of the linearized theory? There 
seem to be at least two ways, summarized in the paragraphs [A] and [B] below, and 
discussed in detail in Section 3. 

[A] The non-perturbative halo that extends to renormalizes the stress-tensor 
of the source, as seen from distances larger than r*. Hence, there is a question of 
matching of the linearized modes obtained for r > r,, to "non-perturbative modes" 
that could potentially be calculated at r < r*. A priori, there is no reason to 
believe that all the linearized solutions can be continued inside the r^, region, and 
yet this is an assumption made implicitly in the works claiming the instability of 
the self-accelerated background. It may well be that certain perturbative solutions 
of the exterior region cannot be smoothly continued inside r^,. An evidence in 
favor of this is a solution of ^2]; in which the interior non-perturbative metric 
interpolates smoothly into a non-perturbative solution outside r*, but not into the 
naive linearized solution. Although the gravitational fields are weak in both interior 
and exterior regions, they have non-analytic dependence on couplings/parameters 
of the theory and cannot be obtained in the linearized approximation. 

[B] There is a known phenomenon of "linearization instability" in theories of 
gravity |H2|, see also [SHI EH ESI EE] • It is said to occur when not all the solutions of 
the linearized approximation can be carried over into the full nonlinear theory. For- 
mally speaking, it arises because the nonlinear equations impose certain constraints 
on the results of the linearized calculations. These constraints cannot be captured 
in the linearized theory itself, and the linearized solutions typically fail to satisfy 
them. The fundamental reason why these constraints arise is this: Gauge fixing 
leaves the residual freedom w.r.t. the gauge transformations by the gauge functions 
that coincide with the Killing vectors of the corresponding background. Then, there 
could exist a charge, that is an integral of a Killing vector contracted with a certain 
current, which is a generator of these transformations. Upon quantization of the 
theory, the physical states should be required to be zero eigenvalue eigenstates of 
the charge operator (the charge is zero in the classical theory by construction). Fur- 
thermore, the requirement that this charge be zero order-by-order in the weak field 
expansion, imposes in the second order perturbation theory an integral constrain on 
the purely first order results. 

The above described phenomenon is rather generic; the only condition it requires 
is that the charge itself be well-defined (e.g., the integral be convergent). This is 
certainly the case if the background space-time has a compact spatial section, then 
the integral in the definition of the charge is w.r.t. this compact 3-space and its 
convergence is guaranteed. In this case, the physical meaning of the additional 
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constraint is simple - a compact 3-space cannot support an isolated "gauge charge" , 
since its flux lines have nowhere to go. 

On the self- accelerated background, similar global constraints can be found in 
its Euclidean continuation. In this case, the analog of a compact 3-space of the GR 
examples is the S"' of the euclideanized worldvolume de Sitter (dS) space. However, 
there is also an additional nonlinear equation on 4D quantities that exists in DGP 
but does not arise in GR. Complementing this equation with the compactness of S^, 
one finds a severe global constraint on linearized perturbations, as we will show in 
Section 3. 

At the moment of this writing, we do not know whether the nonlinear effects 
discussed in paragraphs [A] and [B] improve or worsen the stability, however, the 
following is certain: Any claims on the (in) stability of the self-accelerated solution 
are unjustifiable without studying non-linearities. We emphasize these arguments 
here partially in response to the dissemination of an implicit (and erroneous in our 
opinion) assumption that the question of (in) stability in this model can be addressed 
just by looking at the linearized perturbations. 

The right way to do calculations is different from what's usually done in GR. 
For a given cosmo logical source/overdensity, as a first step, one could calculate the 
perturbations in the domain r ^ r*. There, it is the metric due to the source that 
dominates over the cosmological background. These perturbations are the most rel- 
evant ones, e.g., for the structure formation. One could use the small "graviton mass 
expansion" about the background of the source to see the presence of (in) stability of 
the theory. That the leading order results should be close to those of GR is implied 
by the expansion itself, however, this cannot guarantee the consistency of the results 
(i.e., absence of ghost and other instabilities) in the sub- leading order. These are 
the checks that should be done. 

The above comments also concern attempts in the literature to contrast the 
predictions of the model with observations. In some of those cases, an assumption 
is being made that the values of the cosmological parameters can be extracted 
from the early cosmology using conventional GR, and those results could then be 
used to constrain the late time acceleration in DGP, see, e.g., jTH]. This however 
is not well justified: The structure formation in DGP, although is expected to be 
qualitatively close to that of GR [Hj, is nevertheless influenced by all the non-linear 
effects described above, and will differ from the GR results when it comes to the 
precision tests (see, e.g., comments in 1 1 Tj ) . 



In the second part of the paper we discuss the DGP model with non-conformal 
sources which is artificially truncated at the linearized level. The resulting theory is 
ruled out already by the Solar system data, and presents only an academic interest; 
yet this is the model implicitly adopted in the works which ignore nonlinear effects 
while discussing the question of stability. We show that even in this model there are 
issues that have not been clarified. The main one is that of the boundary conditions 
and poles of the Green's functions from which the content of the physical states of 
the theory can be read off. This involves subtleties which we discuss in Section 4. 
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We do not discuss the quantum gravity loops in the present paper. The loop 
corrections lead to the breakdown of the perturbation theory, like the tree level 
diagrams do. The scale of the breaking is rather low, and one may worry about 
the higher dimensional operators that can be generated at that low scale and spoil 
the predictivity of the theory [201 1^ • The self-accelerated solution itself is right 
at that breaking scale, and any perturbation about it, exceeds that scale. If the 
above approach is adopted, then the analysis showing that on the self-accelerated 
background the bending mode flips the sign of its kinetic term, and, hence, produces 
a ghost [21], cannot be conclusive. These issues, and ways to address them consis- 
tently, were discussed by Dvali [SZl ISH] • In the present paper we adopt the approach 
[T^ in which all the loops are calculated about some external classical background 
field (e.g., early cosmology, galaxies, planets etc). In this case, the loop effects are 
additionally suppressed by the physical scale of the background metric, and can be 
harmless at observable distances. 



2 Perturbations of empty self-accelerated brane 

In this Section we discuss an instructive setup with an empty self-accelerated brane. 
Gravity is truncated to linearized perturbations about the background that is a so- 
lution of the full nonlinear equations. There are no sources placed on such a brane; 
the cosmic acceleration is triggered and maintained by the gravitational sector itself, 
and we look at the evolution of small gravitational perturbations on such a back- 
ground. The obtained results will also be applicable to the theory with conformal 
sources. 

The DGP gravitational action takes the form pP 

5 = ^ y d'x^.Rig,) + ^ j d'xdy^nig) +..., (2) 

where R and TZ are the four-dimensional and five-dimensional Ricci scalars, respec- 
tively, and M^, stands for the gravitational scale of the bulk theory. The analog of 
the graviton mass is rric = 2M^/M|j. The (4 -|- 1) coordinates are x*^ = {x^, y), 
/X = 0, . . . , 3; (^4 and R are the absolute value of the determinant and 4D curvature 
for the 4D components Qi^iy^x^y) of the 5D metric gAB{x,y) 

9f^uix,y = 0) = gAf^A^) ■ (3) 

There is a boundary (a brane) at y = and Z2 symmetry across the boundary is 
imposed. The boundary Gibbons- Hawking term is implied to warrant the correct 
Einstein equations in the bulk. The matter fields, that are omitted in the present 
section, will later be assumed to be localized on a brane. 
Below we consider the metric perturbations of the form: 

ds"^ = {9txv{x, y) + 5g^,y{x, y)) dx^'dx" + 25g^5dx''dy + (1 + Sg^^ix, y)) dy^ , (4) 
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where the background quantities will be denoted by the bar g^u{x, y) = A'^{y)'yfj^^{x). 
We start with the 5D equations of motion. They decompose into {/x5}, {55} and 
{l^jy} components. Let us first look at the {/x5}, {55} equations. They read respec- 
tively 

V^K^, = VuK, (5) 
R = K^- K^.K^'' . (6) 

Here, K = g^^K^y is the trace of the extrinsic curvature tensor defined as follows: 

= ^ {dyg,y - - V.N,) , (7) 

and is a covariant derivative w.r.t. the metric g^^. We introduced the lapse 
scalar field A^, and the shift vector field according to the standard rules: 

g,5 = N, = g^yN- , g,, = + g,, N>^ . (8) 

Equations (0) and © should be satisfied everywhere in the bulk as well as on 
the brane. The {^iv} equation, on the other hand, takes different form depending 
whether it is in the bulk or on the brane 

GifJ =0 for y^O, (9) 

Tfi 

G^.|,=o -^[K,y-g,yK]tl = T,u{x). (10) 
Spectrum of empty self-accelerated brane 

We begin by analyzing small perturbations about an empty (i.e., T^,y = 0) self- 
accelerated brane. We introduce the following self-explanatory notations 

gtiv = (jfiiy + Sgfiiy , = + 5K,y , (11) 

K,y = ng,y, K = r''K,y = 4n n=^, (12) 

where, according to the above definitions: 

K = g^^K^y = K + 6K, 6K = r^SK,, - nr^Sg,, . (13) 

For the self-accelerated solution | = ttLc. The unperturbed part of the off-diagonal 
equation (jSj) is identically satisfied. While the unperturbed parts of the {55} and 
junction equations respectively give: R = 12H^ and R = UHmc, these are consis- 
tent with each other (from now on we will replace \H\ by H with the understanding 
that we should always take a positive value for the latter no matter what is the 
direction of the time arrow). 
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Let us analyze the small fluctuations described by the above system of equations. 
The perturbations of the off-diagonal equation (0) satisfy the following relation 

{6K^, - n6g^,) = V. {g'^SKp^ - ngp'^dgp,) , (14) 

where is a 4D covariant derivative constructed out of g^^. We will return to this 
equations shortly. 

For small perturbations of the {55} equation (jUj) we find 

6R = Gnir^'SK^, - ng'^^dg^,) , (15) 

while from the source-free junction conditions one gets 

5R\o+ = 3m,{r''SK^, - nr''Sg,,)\o+ . (16) 

Since nic = H, equations (fTHjl and (fT?)|) are in contradiction unless the r.h.s. of (fTHjl 
is zero. Moreover, requiring the continuity of the 4D curvature perturbations 6R we 
conclude that the r.h.s. of (fT3j) should also be zero for arbitrary y: 

r''5K^,-nr''5g,, = 0. (17) 

We will come back to this relation below. 

Furthermore, we need to study small perturbations of the bulk {/iz/} equation. 
Since we are dealing with an empty brane (or the one with a conformal stress tensor), 
we can adopt the Gaussian normal coordinates and choose the following "gauge" ^ 

V''Sg,, = 0, r''Sg^. = 0, 6g55 = 6g,, = 0. (18) 

In this gauge the bulk {/xz/} equation takes the form: 

dlSg^u + ^ (^4 - 4i/2) Sg^, = . (19) 

The junction condition reads 

- i (^4 - AH' + Hdy) %,|o+ = . (20) 

Let us discuss the KK modes. We use the following decomposition: 

5g,,{x,y) ^ j h^;::\x)Uy)dm, {D, - 2H')h^;:\x) = m2/i(™)(x) . (21) 

It is useful to turn to a new coordinate z and a function Um 

dz^^^, A^l-'u^ . (22) 



^In general, this is not an acceptable gauge fixing condition for the metric that couples to a 
non-conformal source. That is why one eventually needs to introduce the brane bending mode 



HOI. 
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This choice is very convenient since the bulk equation (fT^ and the junction condi- 
tions (^01) for these variables reduce respectively to Eqs. (jl28|) and (|129|) of Appendix 
A, where the latter equations have been solved (see also |22I)- The spectrum consists 
of a normalizable localized mode of mass ml = SrricH — = 2m^ and a continuum 
of KK modes with masses mj^j^ > 9if^/4. Note that because of the boundary con- 
dition that depends on the KK mass, the modes themselves are not orthogonal. The 
physical modes are their linear superpositions. However, this won't be important 
for our discussions. 

Besides these states, there is a zero- mode solution, m = 0, /o = (see Ap- 
pendix A for details), supported by the above system of equations. This solution 
is not plane-wave normalizable. Its wave-function grows in the bulk as fast as the 
background metric itself. When the warp-factor is scaled out, the zero-mode has 
a constant profile in the y direction, A'^{y)'yfj_i,{x) + A'^{y)h'^^{x), where 'y^u{x) de- 
notes the 4D dS metric and h™{x) is the zero- mode wave-function. We will discuss 
properties of this mode in Sections 4, but we ignore it for the time being. 

Let us concentrate on the lowest massive localized mode mentioned above. The 
mass of this mode is related to the background space-time curvature as follows: 

ml = ^ = 2H\ (23) 

It has been known for some time that 4D massive gravity on a dS background with 
the relation exhibits an additional local symmetry [Hn] 

6h;, = {V^V, + H^^^,)a{x), (24) 

with an arbitrary a. This enables one to gauge-remove one out of the five on- 
shell polarizations of the massive graviton. The resulting graviton has four physical 
polarizations Advantages and disadvantages of this theory are summarized in 
Appendix B. For further convenience, we will denote the Lagrangian of that theory 
by The question that we'd like to discuss is whether the DGP Lagrangian 
on the self-accelerated solution reduces at low energies to [221 123- To answer 
this question we note that the local symmetry (j24j) is present in DGP only for the 
restricted gauge functions. Indeed, both bulk and brane equations of motion are 
invariant, in the absence of the sources (or with conformal sources only), w.r.t. the 
following symmetry transformations 

^ 5g^^ + (V^V^ + H'^i^u)(J , (25) 

where the field a satisfies 

(□ + m"^) a = . (26) 

One can see this by performing straightforward variations of the equations presented 
above and using the following identities: 

(□-4i72)(v^V, + i/V) (scalar) = (V^V, - i/V)(° + 4^^') (scalar) , 

V^(V^V^ + //V^) (scalar) = V^(n + 4//^) (scalar) . (27) 
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In the absence of sources, the symmetry can be used to remove on-shell a poten- 
tially dangerous degree of freedom, and fluctuations can be quantized consistently 
without any ghost or other pathologies, as we'll see in the next subsection. 

Before we turn to these issues, however, we'd like to clarify the following point. 
The difference between the low energy DGP and the model is not entirely clear 
from our discussions above. Indeed, the profile of the localized mode m^, is known, 
it's 4D effective Lagrangian can be calculated by ignoring other modes and integrat- 
ing the 5D action w.r.t. the y coordinate. Doing so, it would naively seem, that 
the low-energy effective Lagrangian on the self-accelerated branch coincides with 
C^. However, as was shown in Ref. j2Sl, the model, in general, does not reduce to 
the theory. This is because of the following: In the calculation of the present 
section the brane bending mode was ignored, it was simply put equal to zero. This 
is consistent for an empty brane (more precisely, as long as T = 0). As a result 
of this, a certain constraint, which is present in the full linearized theory, is not 
captured. To find this constraint, and to see whether it distinguishes DGP from 
one should work with the bending mode manifestly present in the action, as was 
done in [23^]. However, as we will show in the next subsection, the bending mode 
appears as a Lagrange multiplier. As a result, in the absence of sources, both the 
low-energy Lagrangian of DGP and can be consistently quantized to contain four 
on-shell physical polarizations without ghosts or tachyons. 

Effective Lagrangian approach 

Let us start with an illustrative model discussed in The Lagrangian contains 

one gauge field and two scalars, ipi and ip2'- 

C = -\fI - ^(9,y.i + eA,f + |(9,¥.2 + eA,f , (28) 

where are two positive constants. The above Lagrangian is invariant under 

the local symmetry transformations 5ipi = 6ip2 = ^A^ = —d^a/e. These can be 
used to gauge-remove one of the two scalar fields, e.g., one can set ip2 = 0. Following 
[231, consider the case 

/i = f2 , (29) 

which presumably mimics physics of the self-accelerated solution discussed in the 
present paper^. For this particular choice the Lagrangian reduces to [2!^ : 

^ = -\f'^. - Y^d.^^y - /i^^'^mV^i • (30) 
^The authors of |23i attribute this model to R. Rattazzi. 

^The case fi ^ /| is supposed to mimics the physics of a brane on which additional tension is 
introduced. This is more subtle and will not be considered in the present work. 
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By calculating the Hamiltonian of the helicity-0 sector of the above theory pUj). this 
model was claimed to have a ghost j2Sl- 

To look more closely into this issue, let us perform in ipUl) the field redefinition 

B^ = A^ + ^d^ip,, \ = fle^i. (31) 

In terms of the new fields the Lagrangian ()30p takes the form (up to a total deriva- 
tive) 

C = -\Fl + \d^B,, (32) 

where F^j^ is a field-strength for the -B^ field now. This is a Lagrangian of a gauge 
field for which the Lorentz gauge-fixing condition d^B^ = has been introduced 
by means of the Lagrange multiplier A. It is well known how to quantize this 
Lagrangian consistently without the appearance of ghosts or other pathologies. For 
completeness we present here the summary of the path integral quantization, as well 
as the canonical quantization due to Nakanishi and Lautrup (NL) l3Uj . 
In the path integral approach we find 

J [dB][dX]e'^'''''^ = J [dB]S{d''Bf,)e'^'''''^-^^'-\ (33) 

The r.h.s. of the above equation is nothing but the path integral for QED with 
the Lorentz gauge-fixing condition where the Faddeev- Popov determinant (ie)f:|n|, 
which is trivial in this case, is dropped. Needless to say that this path integral 
does not propagate any negative norm states (ghosts) and can be straightforwardly 
completed to a BRST invariant theory (the Faddeev-Popov ghost in this case are 
decoupled) . 

Counting of the on-shell degrees of freedom in the Lagrangian formalism is also 
straightforward; (jH^ imposes the Lorentz condition on the gauge field, removing 
one degree of freedom, yet allows for residual gauge transformations 6B^ = d^a, 
where Da = 0. This can be used to constrain on-shell the second degree of freedom, 
rendering two out of the four modes of B^. 

In the Lagrangian quantization of ()32|) the subsidiary Nakanishi-Lautrup condi- 
tion should be imposed on physical states of the theory: A*^^^|Phys >= (here A*^"-* 
denotes the negative frequency part of the field A)*^. The theory defined this way 
has no perturbative classical or quantum ghost-like instabilities P^ I3Up. Moreover, 

^The NL subsidiary condition is an analog of the Gupta-Bleuler condition in the conventional 
Lagrangian quantization of electrodynamics. 

^We point out that the Lagrangian H32() . at the classical level, could have solutions that are 
not present in conventional electrodynamics AX ■ However, these solutions form different super- 
selection sectors as they imply the presence of some additional physics (e.g., sources at infinity, 
different boundary conditions, etc.). We won't consider such backgrounds here since our goal is to 
study small perturbations on a given background (in particular, on the self-accelerated background 
in the gravity case discussed below). 
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although the canonical Hamiltonian of the theory does contain negative terms, the 
physical Hamiltonian that acts on states of physical Hilbert space, which obey all 
the constraints and NL subsidiary condition, is in fact positive semi-definite. 

Let us now turn to the effective Lagrangian approach to the empty self-accelerated 
brane. The four dimensional effective action for the localized and brane bending 
modes was derived in [23j, which we follow below. The effective Lagrangian takes 
the form 

0/72 

= V) - HA^''{V,V, - 7^.n - 3i/ V)'^ - + (34) 

where A^y = 5g^y{x,y = 0) denotes the tensor mode and ip is the bending mode. 
Furthermore, denotes the Lagrangian of a linearized massive gravity on 

the dS background when the graviton mass and the cosmological constant satisfy a 
special relation ml = 2H^. The explicit expression for C^,{A^i,) is given in p:2j (see 
also Appendix B), but what is important here is that C^{Afj_^) is invariant w.r.t. the 
local transformation, 6Afj_y = (V fj,V u+ H'^1nu)o:{x) , for an arbitrary a{x). As a result 
of this symmetry, does not retain any physical scalar degrees of freedom - 

the only on-shell modes described by are 2 polarizations of a transverse- 

traceless (TT) tensor and two polarizations of a transverse vector [20] However, 
the last two terms of C^s are not invariant under the above local transformation. 
Hence, the question arises whether this reintroduces back the scalar polarization 
that could be a ghost. 

We will show below that the results are very similar to those of a simple gauge- 
field example j2S] considered at the beginning of this subsection. In particular, the 
non-invariant terms in Ces can be recast as a product of a Lagrange multiplier with 
a gauge-fixing condition. Just as in the toy example, small perturbations of the 
resulting theory can consistently be quantized, without any ghosts. To see this, we 
perform a redefinition of the A^,^ field 

V = + (V^V. + H^^^,) uj . (35) 

This change of variables does not affect since C^:{A^^) = C*{B^y). However, 
the last two terms of £efr in (jH^ do get modified. To find the result of this 
let us now introduce the traceless operator P^^ = V^Vjy — |7^,^n, that satisfies 
P'^'^P^^ (scalar) = !□(□ + 4i72) (scalar), then, 

= £*(P^.)-iJP^'^(v^v,-7M.n-3ifV)<^-^<^(n + 4^^')v 



+Huj3H\n + AH^)ip, (36) 



""^^Notice that for ml = 2H^ there exists a unitary massive spin-two representation of the dS 
group with 4 degrees of freedom. 
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and setting uj = -^f, the effective Lagrangian becomes 

£efr = /:*(5^.)-^//(V^V,-7M.n-3i/V)5^'- (37) 

Remarkably, the bending mode (p acts as a Lagrange multipher imposing a certain 
constraints^. This is very similar in spirit to the NL form of the electrodynamics 
Lagrangian (j221)- Hence, our consideration below will be similar too. We start with 
the path integral quantization of this theory 

J [rfi?][rf<^]e^/'^'^^^« = J [di?]5((V^V. - 7M.n - 3H\,)B^ne'^ , (38) 

the r.h.s. of this can be regarded as the path integral of a theory described by 
in which the gauge fixing condition 

(V^V, - 7M.n - 3ifV)^^' = 0, (39) 

is imposed. Moreover, it is straightforward to calculate that the Faddeev-Popov 
determinant for the gauge condition (jH^ takes a form, (iet|n + 4if^|. This is a 
trivial factor (independent of B) which can be restored in the r.h.s. of (jHHj) at the 
expense of an arbitrary overall prefactor. The result reads 

J [dB]S{iV^V, - 7;..n - 3H^'y^,)Bndet\n + 4H^\e'^ '^'^^^ . (40) 

This is nothing but the path integral for the theory with the gauge fixing and 
Faddeev-Popov determinant terms. Therefore, we conclude that the quantization of 
£eff is identical to the quantization of and both of these theories (in the absence 
of sources) propagate 4 physical on-shell polarizations and no ghosts. 

Let us also count degrees of freedom in the Lagrangian formalism. We first recall 
how things work in the theory described by The equations of motion and Bianchi 
identities lead to the condition V^B^^, = V^B (see Appendix B). This equation, 
as well as the Lagrangian C^, itself, are invariant under the local transformations 
5B^y = (V^Vy + H'^'~^^y)a{x), with an arbitrary a. This can be used to gauge 
remove the trace, i.e., set i?^ = 0. As a result, we also get V^B^^ = 0. Hence, 
the B^i, field is transverse-traceless, hence, it has only 5 independent components. 
However, there still remains gauge symmetry under the transformations 6Bf^j^ = 
(V^Vy + H'^'jfj,^) cr(a;), where (□ + 4:H^)a = 0. This allows to eliminate one more 
on-shell degree of freedom. Hence, 5^,^ field in the Lagrangian £^ describes only 4 
physical polarizations 

Let us now see how the counting works for C^s- Variation w.r.t. (p yields (jSHI). 
Furthermore, it is straightforward to check that the Bianchi identities for the {yuz^} 

"Notice that (V^V^. - j^,^a - 3H^-ff,^)Bt''' ~ SR. 

^^The Lagrangian C^,, however, does not permit the Bf^i, field to be coupled to a non-conformal 
source, T ^ (see e.g.. Appendix B). This is very similar to how the conventional electrodynamics 
does not permit the photon to be coupled to a non-conserved current. 
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equation lead to V-B^^ = V yB. The divergence of the latter combined with 
leads to 5 = 0, and, as a consequence, V^B^^ = too. Moreover, the Lagrangian 
Ces, as well as all the above constraints, are still invariant under the transformations 
6B^y = (V ^Vu+H'^'y^u) cr{x), where (□+4if^)(T = 0. These can be used to eliminate 
one more degree of freedom of the Bf^^y field, reducing the total number of its physical 
on-shell states down to 4. Therefore, as far as the B^^, field is concerned, there is 
no difference between the description in terms of or C^s- Let us now look at the 
(p field. This field is a Lagrange multiplier in (j37p . Moreover, the trace of the {/xi^} 
equation plus the Bianchi identities lead to (□ + 4H'^)(f = 0. This is reminiscent 
of the condition DA = on the Lagrange multiplier emerging in ()32j). As in that 
case, the value of the Lagrange multiplier (p can be set to zero and this condition 
should be possible to maintained in the full quantum theory by imposing the NL 
like constraint on the physical states y9'^~)|Phys) = 0. This must be sufficient to 
guarantee that the small perturbations described by the Lagrangian Ces{,B) yield 
no ghosts, as it was shown in the path integral quantization. Moreover, the physical 
degrees of freedom are identical to those described by 

Summarizing, we showed that small perturbations on the empty self-accelerated 
brane can be quantized consistently and do not exhibit any ghost instability. These 
results also apply to the brane endowed with conformal sources. This is because all 
the arguments and counting of the degrees of freedom go through in the presence of 
conformal sources, since the symmetry ()25 |) . (Pm) is preserved in that case^^. 

The above presented counting of the degrees of freedom in Ces changes once 
sources with T 7^ are introduced. Details are given in Section 4; here we just note 
that although the relation V^-B^z^ = V^B remains intact, one looses the symmetry 
fl25|) ■ (pUj) that was essential to remove one on-shell degree of freedom. The physical 
metric that couples to the stress-tensor is 

= Sg^,, + ^^(V^V, + H\,Mx) . (41) 

Moreover, {n + 4H^)ip ~ T, and a nonzero-trace source could excite a ghost and be- 
come unstable. However, such a source also leads to the breakdown of the linearized 
approach adopted above. We turn to the discussions of these issues in Section 3. 



3 Perturbations in the presence of sources 

As we have already discussed in Sections 1 and 2, non-conformal sources could excite 
ghost in the linearized reduction of the theory (see, however, discussions in Section 
4). If this were all, the sources would start to radiate ghosts until they'd disappear 

""^^In this paper we do not discuss stability of the background w.r.t. possible classical, and semi- 
classical non-perturbative effects, such as the decay of the vacuum etc. If these instabilities are 
present, we'd expect them to set in very slowly, without affecting the observations; moreover, the 
full nonlinear theory should be used to study their consequences. 
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(if the sources carry no conserved charges), or would accrete some negative "gradient 
energy" to produce an effective source that by itself is conformal. However, in all 
these considerations it is important to remember that non-conformal sources lead 
to the breakdown of the linearized approximation; they surround themselves with a 
nonlinear halo of a variable curvature that locally dominates over the self-accelerated 
background. The naive linearized solutions, that could contain instabilities, are not 
guaranteed to penetrate the halo. We will present an evidence that the linearized 
solutions are indeed spurious in this case. Dynamically, the halo screens the bare 
mass of the source, and because of this, the process of its formation is energetically 
favorable. 

Complete set of nonlinear equations 

To discuss these issues it is useful to start with the DGP action in the ADM for- 
malism gS] (see, e.g., 

S=^l d'xdy^g{R5{y) + '^N [R + - K,,K^^)) . (42) 

Equations of motion of the theory are obtained by varying the action (jl^ w.r.t. 
g^u, N and A^^. Here we reiterate a subset of two equations, the junction condition 
across the brane, and the {55} equation that can be obtained by varying the action 
w.r.t. N . The former reads as follows: 

G^, - m,{K^, - g^,K) = T^jMl, , (43) 

where G^j^u is the 4D Einstein tensor of the induced metric g^j^y and T^^ is the matter 
stress tensor. The {55} equation takes the form 

R=K^-Kl. (44) 

Note that (^Hj) is valid only at ?/ = 0^ while should be fulfilled for arbitrary 
y. The terms with the extrinsic curvature contain derivatives w.r.t. the extra 
coordinate as well as the {55} and {/i5} components of the metric. By expressing 
from K and K^j^y in terms of Rfj^y and T^j^^ we obtain the following expression: 

,nlR 3Rl^) + J, (2R,.T^ - Irt) - J, (t^I - ir=) ,(45) 

The remaining equations correspond to the {/i5} and bulk {fiiy} equations. The 
{/i5} equation for arbitrary y (which we copy here again for convenience) reads as 
follows: 

V^K = V'K^, . (46) 

The covariant derivative in the above equation is the one for g^^i,. At y = Eq. 

is trivially satisfied due to (jlHjl . For y ^ 0, (jlUj) sets the relation between A^"^, A^ 
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and Qf^^. Hence, (HUj) gives a relation between these quantities for both y = and 

One can use the bulk {55} equation (j44p to determine N in terms of A*"^ and 
Qfj^iy, and then use (pUj) in order to express A^^ in terms of g^u- If so, there must 
exist one more equation which should allow to determine the bulk g^^, itself. This 
is the bulk {/xi^} equation, to which we turn now. The latter can be written in a 
few different ways. For the case at hand, the formalism by Shiromizu, Maeda and 
Sasaki j^Hl is most convenient. In this approach, the bulk {/xi^} equation and the 
junction condition can be combined to yield a {/iz^} equation at y ^ 0^. This gives 
a "projection" of the bulk {/ii^} equation onto the brane. Since the bulk itself is 
empty in our case, this is the most restrictive form one can work with. The equation 
reads as follows gHl 113 UHl Hi : 

Gnu = — ^Bav — Enu , (47) 

where 

B^u = —G^aG'^ + '^G^G^u + -g^uGaisG'^^ — q9i-iu{G'^)'^ , 

G,u = G^u-T^JMI,, (48) 

and where all quantities are taken at ?/ = 0^. G^j.^ denotes the 4D Einstein tensor of 
the metric g^y, E^^ denotes the electric part of the bulk Weyl tensor, projected onto 
the brane. An important property of the Weyl tensor is that it is invariant under 
conformal transformations. Moreover, E^j^ is traceless. 



Perturbations in nonperturbative domain 

One way in which the presence of instability was argued in the literature is 
by looking at the scalar bending mode and observing that its kinetic term flips the 
sign due to curvature of the self-accelerated background. However, the scalar mode 
carries no gauge- independent physical information^^. Because of this, no conclusion 
can be obtained just by looking at this mode; as we have seen in Section 2, the 
mixing with the tensor should be taken into account. 

To emphasize once again how misleading the scalar mode could be, we briefly 
recall below a well-known example. Consider conventional 4D GR coupled to matter 
on a fiat space. No instabilities are present in this theory. Let us perform now 
conformal transformation of the metric by the conformal factor (I+tt) . The resulting 
Lagrangian contains the tensor field, kinetic mixing term of the tensor to vr, and, 
most importantly, a wrong-sign (ghost-like) kinetic term for the scalar vr. Moreover, 
this scalar mode does couple to matter in the new frame. If one simply ignores 

^"^This is true even in the so called decoupling limit where as was shown in 14, the scalar mode 
does not decouple in DGP from the tensor modes. 
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the tensor mode (by putting it equal to zero, or by considering fields on which the 
conformally rescaled R is zero), one would naively conclude that tt is a ghost that 
could be emitted by the source. However, this would be an erroneous conclusion, 
since we started with the theory with no ghost in the first place. The error is made 
by neglecting the mixing of vr with the tensor, and ignoring the fact that vr by itself 
carries no physical information. 

The above example is instructive, and does certainly demonstrate the danger of 
drawing conclusions in gravity (as well as in gauge theories) based on the helicity-0 
sector alone. However, the situation is not as simple in the DGP model. This is be- 
cause there are nonlinear mixing terms between the scalar and tensor modes which, 
on classical backgrounds, give rise to quadratic mixing terms between them. Because 
of this, as we will argue below, what determines the mixing at scales r ^ r*, is the 
local background due to the source, and not the global cosmo logical backgrounds^. 

Let us consider small perturbations around a classical solution (e.g., such as the 
nonperturbative solutions of ^H]) where there is significant curvature for distances 
smaller than a certain distance r^:. Taking perturbations of the following form: 

RpLv = R% + ^Rtiv , (49) 
we obtain from the first order expansion of (j45p 

mlSR = lR''SR-'^R>R''' + 4l,^tT^^'^-^R''T. (50) 

Notice the crucial role played by the non-vanishing curvature of the classical solu- 
tions in determining the properties of the perturbations. 

Let us first discuss the conventional branch. In this case, for distances much 
greater than the Ricci tensor is negligible. This yields 6R = (we consider 4D 
sources localized to a region smaller than r^). While in the regime of distances 
smaller than r^, we get 

ml6R = Iw'^HR-2R^^^5R^^ + ^W^^T^'^ ^ (51) 

where the superscript "sch" refers to the nonzero curvature solution in the conven- 
tional branch. 

On the self-accelerated branch, in the large distance regime where the curvature is 
de Sitter like (with Hubble parameter equal to rric) Eq. (jSOI) reduces to 5R = —26R, 

^^It is worth mentioning that the dynamics of smaU perturbations of the bending mode tt of Refs. 
|2l)[l21j . in the nonperturbative regime, is completely dominated by the local background. That is 
to say, the fluctuations of the tt field are well above the strong coupling scale of the theory, the later 
being determined via the graviton scattering amplitude on Minkowski background [20]. However, 
the "effective" strong coupling scale of the theory, even in the case of the graviton scattering on an 
empty space, is determined by the local non-perturbative background due to gravitons themselves, 
or due to any other supercritical sources which are present in a given setup. This effective scale is 
much higher than the strong coupling scale of [201 (see, discussions in the first reference of |13|). 
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which amounts to 6R = as in the conventional branch. In the short distance 
regime the curvature flips the sign with respect to the conventional branch^^: 



- -^p • (52) 

Therefore, yields, in this case 

ml6R = -Iw^HR + 2T^^5R^^ - ^T^^T^^^ + (53) 

with only the sign differences of the r.h.s. as compared to ()51|). 

At scales below r,, the local curvature R^'^^ overcomes m^, and becomes larger 
and larger as one approaches the source (still remaining much smaller that the 
magnitude of the Riemann tensor components of the conventional 4D Schwarzschild 
metric ^H])- Because of this, the expressions on the l.h.s. of (j^Tj) and (j^Hj) can be 
neglected as compared to the terms on the r.h.s. In this approximation, the above 
two equations coincide up to an overall sign. 

Very similar calculations can be performed with a more general equation (j47|) . 
Since the conclusions of those studies are analogous to what we found above, and, 
moreover, the above presented expressions can be obtained by taking trace of (|17j), 
we won't present those results here. 

Therefore, the dynamics of small perturbations at scales r ^ r^, is completely 
determined by the local metric due to the source, and not by the cosmological 
background. The linearized perturbations, which by themselves are valid at r > r^,, 
are not guaranteed to match upon the solutions inside the r ■^r^: domain. Physically, 
this is because the halo of a variable curvature that extends to scales r ~ r^,, screens 
the stress-tensor of the source [T^ , making it hard for the naive linearized solutions 
to match it smoothly. For instance, in the static solution of ^HIj the modes that 
at r > match smoothly the interior solution, are not the modes of the naive 
linearized approximation. The latter are spurious. 

Global constraints on self-accelerated background 

In this subsection we consider Euclidean continuation of the self-accelerated back- 
ground and discuss the global constraints imposed on the linearized perturbations 
in a theory with, and without, sources. 

We start with the constraints similar to the ones emerging in 4D GR on a dS 
space-time with a closed spatial section. This part of our presentation is a straight- 
forward generalization to one more dimension of the approach of Ref. Later 
in this subsection we will discuss the constraints that are more specific to the DGP 
model. However, we should point out that, in terms of the Lorentzian signature 
space-time, the results of the present section are specific to a foliation of the dS 

^^This statement is a direct consequence of the structure of the Israel junction conditions in 
DGP, and is independent of particular details of the solutions at hand. 
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space-time that covers the whole dS hyperboloid. We do not expect similar results 
to arise for flat and open dS backgrounds in the absence of non-conformal sources. 
On the other hand, if the non-conformal sources were present, such global con- 
straints could formally arise for different reasons. However, since in the latter case 
the linearization breaks down, these constraints should have no importance for the 
self-accelerated solution. 

Consider the euclideanized dS space, S'^. Suppose there exists some symmetric 
tensor n^'^ which is covariantly conserved 

V.TT^'^ = , (54) 

with Vj, being the covariant derivative compatible with the metric Qf^,^. Then, the 
integral below is a total derivative, and, assuming that fields are single- valued, it 
vanishes on the S"^ 

Q^= [ d^x^{£^g^,) TT^'^ = . (55) 
The Lie derivative £^ w.r.t. the vector field is defined as follows 

= rv,v+ vv.r +t.aV^r, (56) 

where t^^, is an arbitrary smooth tensor and the second line follows from the first one 
given that is a covariant derivative of an arbitrary metric whose affine connection 
is torsion- free {g^v is one example, but not the only one). 

Let us split g^j_i, and n^'^ in their background and perturbative parts 

g,, = + h^, , TT^'^ = TT^'^ + p'^^ . (57) 

Suppose that is a Killing vector of the background space-time, that is 

£a^^u = . (58) 

Moreover, below we will consider tt'^'^s for which 

£^n>"' = . (59) 

Using these properties we expand the expression for up to the second order 

in small perturbations. Ignoring the surface terms on the S'^, we get 

^ 1 d'x^ {£^h^^l) = . (60) 

Here h^^^ and p^^^ are the solutions of the first order equations. Expression (jUUI) is 
an additional constraint imposed on these perturbations. 
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Let us now turn to the expression for ir'^'^ that are specific to the model at hand. 
There are at least three tensors that satisfy the properties of n'^'^ described above. 
These are: 

Kt'^-gi^^K, (61) 

Only two of these three are independent as they are related to each other by the 
junction condition G^"" - rndK^"' - gf^'K) = T^^/M^^. The global constraint ^ 
for these tensors reads respectfully: 

Qf = j d^x^{£^h^^l){K^'^^^ -{g^^Ki''^) = Q, (62) 
Qc = I d'x^ {£ch(;}) G^'^^'^ = , (63) 
Qc = Jd'x^{£^hl^:i)T^'^>'-' = 0. (64) 

Any one of the above constraints follows from the other two. Let us focus on the 
last one (j64p . In an expanded form it reads: 

J d'x^ (C"VW/i« + /^Svi^^C + /^ti^vWc) T(i)^'^ = , (65) 

where V^a^ denotes the covariant derivative compatible with the background metric 
7^,y. The significance of this constraint is the following, one would calculate h^^^ in 
terms of the vi''^'* and an arbitrary source T*^^^ (as it is done in the next section). 
The question whether the obtained /i*-^-* can or cannot be continued into the full 
nonlinear theory would have been decided based on whether or not it could satisfy 
the constraint (j65p . The latter, in general, is hard to be fulfilled automatically, and 
it would impose a constraint on the acceptable sources T*^^^ However, in the case of 
the self- accelerated solution with the sources the situation is slightly different. Local 
curvature within an arbitrary overdensity above the background is bigger than m^. 
As a result, the perturbative expansion should break down for all these sources. 
Because of this, the derivation of (|65|l presented above fails. The effects of the 
nonperturbative domain lead to screening of the mass of the source jl4j . 

In the remainder of this subsection, we would like to discuss the global constraints 
on perturbations of the self-accelerated background without sources. These are more 
specific to the DGP model, as they are based on a nonlinear trace equation that is 
very different from that of GR. However, as we mentioned at the beginning of this 
subsection, these constraints are specific to the foliation of the dS space-time that 
covers the whole dS hyperboloid. 

To get to the point as quickly as possible we again turn to the euclideanized 
version of the background dS space with the topology of a four-sphere S"^, and look 
at one of the exact equations of motion of the DGP model 

?,mlR = - SRl^ . (66) 
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This equation is valid at the position of the brane in an arbitrary coordinate system. 
We will be expanding this equation about the background up to the second order 
in perturbations (which corresponds to the expansion up to the cubic order in the 
action). Hence, 

g,u = lf.u + Sg^., g,, = Y' -dg'^' + 5g^''5gl + ..., (67) 

where the first equation in (jU7j) is just the definition of 5gij,y. Furthermore, we use 
the following standard notations for the expansion of the Ricci tensor: 

R>.. = 3H'j^, + + ^ (68) 

where R^^i! and R^^u denote respectively the linear and quadratic in perturbations 
terms (for convenience, the explicit expressions for these quantities are given in 
Appendix C, see also [SO])- It is also convenient to adopt the following notations: 

%. = /^!.'j + /^S + -, (69) 

where h^fj^} denotes a solution of the first order equation, while /i^fj is a solution of 
the second order equations. The decomposition of (jU^ is useful for keeping track of 
terms in the corresponding orders. For instance, the term R^^uih^^^) is a first order 
quantity, while the terms R^^}{h^'^^), R^ii}{h^^^) and h^^^^''' R^^^h^'^^) are all the second 
order quantities. On the other hand, the terms h^'^''^" R^i^'uih'^^''), h'^^^^'^ R^^,!{h^^^) ap- 
pear in the third order, while R^^}{h^'^^) in the fourth order. We will be ignoring all 
the terms beyond the second order. 

To proceed further we need to specify the gauge conditions at each order. In the 
first order we choose the gauge 

I'^h^^l = . (70) 

The first order source free equations of motion also allow to set (although the latter 
is not necessary for our derivations) 

V^/iW = . (71) 

In the second order, however, we impose only the tracelessness condition ^^: 

T^'^/^g) = . (72) 

-•^^This does not fix the gauge freedom in the second order completely, but is enough for our 
purposes, i.e., our considerations apply irrespective of an additional gauge fixing condition that 
one may need to introduce to fix the residual gauge freedom of the second order perturbations, as 
long as that conditions is consistent itself. Notice that we do not impose the covariant transversality 
condition in the second order, as this would have over-constrained the system. 
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Now we are ready to proceed to the expansion of ()66|) . Since rric = i7and7^^i?ii^(/i«) = 
, we find in the second order in perturbations: 

R = g^^R^, = 12H^ + Y'R^^Kh'^'^) + T^'^^lH^^'^) 

-/^«^'^i?«(/^^'^) + 3i/'(/^il))^ (73) 

The terms on the r.h.s. of ()66|) are expanded to the second order as follows: 

-24if2/,(iWi?a) (/,(!)) + nH\h^^lf , (74) 

and 

Rl, ^ R.^Ro^pg^^g"^ = 36H' + 6i/%^^<)(/i(2)) + GH'^'^R^^Xh^'^) 

-12H'h^'^>^'^RW(h^^)) + 21H\h^^lf + {R^iKh^'^)f . (75) 

Substituting the expressions (f7H|) . (f7i|) and (f7H|) into the equation (jUU]) we get 

-18i/^(/.W)^-3(i?«(/^«))^ = 0. (76) 

All the terms in the above equation are of the second order. This equation should 
be used to determine the second order metric perturbation /ijfj in terms of the 
squares of the first order perturbation h^^} and its derivatives. However, there is 
more to it: Only the first term on the l.h.s. of (fTUI) contains the second order metric 
perturbation the rest contains only the first order metric /i^/j. Moreover, since 
is traceless, then 7'^'^i?|/j(/i*^2)) = V^V^/i|f2. As a result, the first term on the 
l.h.s. of (fTUj) . which also happens to be the only term containing /i/jj, is a total 
derivative. Therefore, if we integrate both sides of (fTHjl w.r.t. the background space, 
which in our case in S"^, the term containing will drop out as a total derivative. 
The resulting integral equations reads: 

/ d'x^ {H'r^R^^^ih^'^) + 5if^/^(^)-^<)(/.«) - 6H\h^^}r - (i?«(/^«))^) 

= 0.(77) 

This equation contains only the first order perturbations h^f^u, and represents an 
additional constraint on these perturbations imposed by the second order equations. 
This constraint cannot be captured in the linearized theory. Although we derived 
(fTBj) using a particular gauge, it represents an integral of a gauge-invariant equation 

^®In conventional GR on 5^ Eq. H73() can be used to derive a global constraint. 
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Can the first order perturbations that are solutions of the linearized equations 
satisfy the above constraint? To answer this question we use the explicit expressions 
for the quantities in (f77|) (see, e.g., Appendix C), which after integration by parts 
yields: 

/ d'x^ {iUh^^f + hH\W^h'~^f) = . (78) 

The integrand on the l.h.s. of (j78p is a positive semi-definite quantity in the Eu- 
clidean space that we are dealing with. Hence, the global constraint (fTSj) implies 
that nh^i^u = 0, and Va/il/j = 0. These can only be satisfied by trivial solutions! 

Summarizing, the nonlinear equations impose very severe constraints on eu- 
clideanized linear perturbations. If fact, on a brane with no matter sources consid- 
ered above, no non-trivial perturbations appear to survive these constraints. Then 
what is physics of perturbations on such a background? The answer to this ques- 
tion, as we discussed before, lies in non-perturbative dynamics. It so happens that 
realistic sources, whether they are due to nonlinear self-interactions of gravity itself, 
or external matter, give rise to the breakdown of the naive linearized approximation 
used in this section. This breakdown obfuscates the meaning of the conventional 
linearization approximation and constraints discussed above. On the other hand, it 
also suggests that the required calculations should be done using different approxi- 
mations, e.g., the one we briefly discussed in Section 1. 

4 Truncated theory with sources 

In this Section we consider the self-accelerated background on which the theory 
is truncated to the linearized level. As we have discussed before, such a theory 
is ruled out already by the Solar system data because of the vDVZ discontinuity. 
Yet, this is the model based on which the conclusions of the instability of the self- 
accelerated background have been reached in the literature. As we argued in the 
previous Sections, the only model that makes physical sense to discuss is the one in 
which non-linearities are retained, whereas the truncated model has only a marginal 
interest. Nevertheless, it can be used to discuss certain theoretical questions, such 
as its relation to the model of massive gravity with the enhanced symmetry [5^ . 
However, even in this case, there are subtle issues concerning the analytic structure 
and boundary conditions of the Green's functions, which have not be clarified. This 
issues are important for determining the linearized spectrum of the theory, as we 
will discuss in detail in the present section. 

The square root and boundary conditions 

To get started, let us pick a coordinate system in which the brane is straight and 
located at the y = point in the fifth dimension. Let us look at the trace of the 4D 
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Einstein equation at the brane without matter or radiation 



R = SttIcK . 



(79) 



Both sides of (f?^ are calculated at y = 0'^. The key point is that K can take 
negative or positive values depending on the boundary conditions in the bulk. This 
could be seen, e.g., in a gauge in which K^^, oc dyg^^{x,y), where g^^{x,y) denotes 
the {yuz^} components of the 5D metric. If the metric falls off the brane in the y 
direction, then the r.h.s. of ()79|) could be negative. We call this the (— ) branch. 
Alternatively, if the metric grows off the brane in the y direction, the r.h.s. of (f?^ 
could be positive. We call this the (+) branch. Clearly, both the (— ) and (+) 
branches admit Minkowski space-time as a solution of (|79|). However, the key point 
is that the (+) branch, in addition, admits the self-accelerated solution. 

To get more insight into this set of issues, let us rewrite (f?^ as a Friedmann 
equation for a spatially-flat background It is instructive [51J to discuss the 
properties of the (— ) and (+) branches in terms of this equation. For the (— ) 
branch it takes the form: 



where H denotes the Hubble parameter as seen by a 4D brane observer. The so- 
lution H = corresponds to a Minkowski background of the (— ) branch. Small 
perturbations about this space are stable [T]. 

For the (+) branch, on the other hand, the Friedman equation reads: 



There is a set of two solutions. The solution (A+) with H = 0, and (-B+) with 
\H\ = rric- The [A^) solution is just a Minkowski background of the (-I-) branch. The 
solution (-B+), that is called the self-accelerated solution [2], is the most interesting 
one, as it describes an accelerated expansion of the Universe. 

In addition to the freedom in choosing boundary conditions in the bulk there is 
also a freedom in setting a direction of the time arrow t on the brane. This freedom 
is reflected in the above equations as a symmetry H — * —H. This suggests that 
a solution with a positive (negative) H could be expanding (contracting) for one 
choice of the direction of t and contracting (expanding) for the opposite choice. 

Some of the above described properties of the (— ) and (+) branches can also be 
seen by analyzing the perturbative amplitudes on a flat space. For them, the opera- 
tor dy in the expression for K^^, is replaced by the square root of the d'Alembertian 
±V— A choice of the bulk boundary conditions is in one-to-one correspondence 
with the choice of the sign of i^^— ^4. Furthermore, the choice of this sign accom- 
panied with the choice of the direction of the time arrow has to do with the causality 
and unitarity of the amplitudes; it determines whether there could or could not be 
tachyon and ghost instabilities in the linearized theory. Since this is an important 



(80) 



- mc\H\ = 0. 



(81) 
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issue, and has direct relevance to our calculations, we will reiterate and sharpen 
below some of the discussions of Ref . [HI] . 

For the (— ) branch, the one-particle exchange amplitude for two very weak 
sources on the brane (in the momentum space w.r.t. the four worldvolume co- 
ordinates) takes the form ^J: 

T-^^/gCxp ( - ^ ^ 

M^~\p, y) = — — ^ ^ , (82) 



where we denote the Euclidean four-momentum squared by = —pl+Pi+pi+pl = 
pi +P1+P2 +P3, and 

Ti% = 8 TT (t'^^ - 1 T° ■ . (83) 

The amplitude (jH^ was constructed by imposing the decreasing boundary condi- 
tions in the y direction for positive Euclidean momenta. In the Lorentzian signature 
this corresponds to an oscillatory dependence on the y coordinate with a plus sign 
of the exponent ~ exp{+i^y —p'^\y\) . In the physical domain, —p'^ > 0, this 

corresponds to a retarded amplitude for a given choice of the time direction. Fur- 
thermore, there is a branch-cut on the complex plane of — on the positive real 
semi-axes due to the square root in (jH^ . and the pole at a/p^ = —rric is on a sec- 
ond nonphysical Riemann sheet (see, ^] for detailed discussions.) This pole has a 
positive residue, while the residue of the pole at = is zero. Thus, it describes a 
metastable graviton without ghost or tachyon pathologies. 

For the (+) branch, on the other hand, the amplitude is determined by choosing 
the other sign of the square root. As a result, one obtains a Euclidean amplitude 
that grows with y 



T^,^ exp ( 

M^^\p, y) = — ^ . (84) 



— rric \/ p^ 

In the Lorentzian signature this corresponds to an oscillatory solution with a nega- 
tive exponent MS^^ ~ exp(— za/— p^l^/l), and describes an advanced amplitude. The 
expression ()84|) differs from (jH^ also by the position of the pole in the denomina- 
tor. The pole at = rric is located on a first Riemann sheet on a negative real 
semi-axes of — p^, hence this pole is tachyonic. This suggests that the space-time 
flat solution of the (+) branch is unstable (this instability could be stronger than 
tachyonic)! However, the situation can be reversed by replacing the time coordinate 
t, by —t on the (-I-) branch. Only in this case the Minkowski solution of the (-I-) 
branch has a consistent physical interpretation. After this replacement, ()84|) will 
describe a retarded amplitude for the new choice of the time direction. Moreover, 
the tachyonic pole becomes a resonance type pole. This is simply because (jHH) turns 
identically into 
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On the self-accelerated background of the (+) branch one can calculate two 
amplitudes with the behavior similar to (|82j) and (|84p . as we will show in this 
Section. However, unlike the flat-space example discussed above, the replacement 
t —>■ —t does not transform one amplitude into the other one. 

Summarizing, there are two choices to be made, one for the direction of y and 
another for t. On a Minkowski background these two choices are degenerate. How- 
ever, on a curved background the degeneracy is hfted and, in general, one should 
study each particular case separately, as it will be done in this Section. 

Spectrum in the presence of a source 

Let us introduce a source Tab{x, u) = ^A^B'^fj.u{x)6{y). In the gauge that we used in 
Section 2, the source would bend the brane [lU]- Hence, we define a new coordinate 
system {x, y) in which the brane is straight, and impose the following gauge fixing 
conditions 

V''^. = V.r^dg^p , 5^55 = S~g^5 = . (85) 

The new and old coordinate systems, and fields, are related by 5D gauge transfor- 
mations ()41|) . In terms of the new variables the metric is written as 

ds^ = {gt,u + dx^dx" + df . (86) 

Let us look at the equations of motion. The off-diagonal equation (0) written in the 
new coordinate system is identically satisfied. On the other hand, the {55} equation 
(jni) takes the form: 

- ^n'r''5~9,u = Qn (g^'SK^, - nr^Sg,,^ , (87) 

which can be solved to obtain 

r''SK,, = -r''S~g,.. (88) 

The latter is consistent with (j4ip . 

The junction condition, taken at 0+, reads as follows: 

- ^ (04%. - Vf,VJ~g2) + 2H'6~g^, - ^^^Jg^ - H [dK^, - i^^i) = T,, -(89) 

The first four terms on the l.h.s. of ()89p is what appears in a 4D massive gravity at 
the point of enhanced symmetry (see Section 2 and Appendix B). Here, however, we 
get additional terms due to the extrinsic curvature (the last two terms on the l.h.s.). 
This makes the self-accelerated background of the present model different from the 
massive 4D gravity at the special point [22] ■ Because of this, in the present model 
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a non-conformal matter stress-tensor is allowed. Indeed, the trace of the junction 
equation takes the form: 

3m,r''SK^u=T. (90) 

(The l.h.s. of this equations would have been zero for a massive 4D gravity at the 
special point (j^ . see Appendix B). 

One way to solve for the induced metric 6g^i, is, to use (jlH), go to the variables 
without the tilde sign, and solve first for them. For 6g^u we find the following 
equation (valid at 0^) 

- i (04 - AH' + Hdy) 5g,, = T,. - ^^,^T + U,,^^^-L__ , (91) 

where U^i, = (V^Vu + H'^l^j.u), and we used SK^y = 5K^^ + Il^u^- Here we note 
that the last term on the r.h.s. of (|91|) is somewhat alarming. It already has a pole 
(all operators are acting from left to right). As a result, the expression for 5g^j_y may 
end up having a double pole. We will see whether this potential problem can be 
avoided in the expression for the induced metric 5g^v 
The trace of the junction condition (jUUI) determines ip: 

We pause here. The pole in the above expression for Lp is tachyonic, it corresponds to 
a scalar of a negative mass square = —AH'. One could conclude that the brane 
bending mode ip is at least as bad as a tachyon. However, (p is gauge dependent, 
and, in order to reach a conclusion one needs to calculate a gauge invariant physical 
amplitude. 

Finally, for future reference we write the expression that encompasses both the 
bulk and junction equations for 6gfj,^: 

[ndl + + 25(Z/)} (^4 - AH')^ 5g,, = -A (^T^, - ^^^,T + HU,,^^ S{y) .(93) 

We will use the above equation for the determination of the induced metric in the 
next subsection. The reminder of this subsection discusses a certain subtlety that 
emerges in the above equation. A reader who is not interested in these details could 
skip to the next subsection, where the actual calculations are performed. 

This concerns the last term in the parenthesis on the r.h.s. of (jU^ (we will call 
it A, 



A,. = (V,V. + H'^,,) ^^^^f ■ (94) 

Looking at (jMj) one may conclude that it has a tachyonic pole at ^4 + AH' = 0. 
However, this would be a misleading conclusion. To see this we use the identity 

(□4 - AH'){V^V, + i/^T/..) (scalar) = (V^V, - H'^^,)iD^ + AH') (scalar) , (95) 
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which allows us to rewrite A^j, as: 



A.. = □ ( V,V. - i/^,.) I . (96) 



As it is clear from (|96|). the expression for A^^, seems now to have a pole at ^4 — 
AH'^ = 0, for the operators that multiply it from the left. We will carefully deal with 
such subtleties in the next subsection. 



Poles and residues 

In order to find the physical spectrum of the theory one should look at the poles 
and residues of the induced metric: 

2A 

Sg,,u{x, y) = 5g^y{x, y) + —Uf.u'P ■ (97) 

Substituting (jHZI) into (jHH), and using the identity we find the equation 

- ^ (04 - AH' + Hdy) 5~g,, = T,, - ^7m.T - (V^V. - H^^.J) ^ . (98) 

As before, this equation is written for y = 0~^. Note, that the troublesome double- 
pole structure that was present in has remarkably canceled out in the above 
equation. 

In order to find the expression for 5gfj_i, that would capture both the bulk and 
junction equations, we use (jHTj) in (jHS))- The resulting equations for Sg^u reads 

+ { ^ + ^^(y) } (°4 - 4i/2)^ S~g^^ = 

-4 (t,, - ^^,,T - (V^V. - H'^,,) 6iy) + j (V^V. - if^^.) ^ . (99) 

Here we would like to emphasize that the last term on the r.h.s. of ()99p does not 
vanish outside of the brane. Therefore, even a source that is localized on the brane, 
gives rise to an effective source that is delocalized all over the bulk. It is the effective 
source that excites the physical metric 5g^y. Moreover, this effective source grows 
in the bulk (this is because in physical amplitudes the last term on the r.h.s. of 
fj99|) should be multiplied by an additional factor of A' arising from the norm of 
the inner product). This rather unusual behavior is a consequence of the fact that 
in the coordinate system in which the brane is static and straight the background 
itself grows in the bulk as A'. 

Let us now solve (j^HI)- It yields two independent solutions 

_ 1 .a±f0W) 1 I ni^) ( T T^_n-^ 

^og^u± ^(4)^ O^^) - H'^a±{0^^))\ V ^ J 3 

^ ^ (100) 



3i/2 (9(4) - 
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where 0^^\ H^^ and a± are respectively defined by 

(101) 
(102) 

(103) 

It is important to point out that the r.h.s. of ()100|) has potentially two poles: the 
first one for the zero eigenvalues of O^^^ and the second for zero eigenvalues of 
(9(4) — H^a± (O*^^^). These poles and their residues will be discussed in detail below. 
Equations ()100|) can be rewritten as 

= ^±((^^'))[tJJ + 1(p,.-^-^p,.)t 

^ + -i-7,.T , (104) 



K±{z) = (105) 



3/720(4) 12H2 
where P^, = V^V. - \j,,n, {i.e., H;, = P^, - li^.O^"^^), Q(^) = -04 - ^H^ 

z — H'^a±{z) 

and the decomposition of a covariantly conserved symmetric tensor in terms of its 
transverse-traceless, spin-0 longitudinal, and pure trace parts was used, namely, 

T,. = TjJ + ^7m.T + Ip,u^T. (106) 

Noticing that 

p^,(n4 + 8H^)ip = n^p^^if , (107) 

for any scalar (which implies that P^yf{Q,^^^)Lp = f{0^^^)P^yLp for any smooth 
function /) the term in the parenthesis in ()104j) vanishes. 

Next, let us introduce the Lichnerowicz operator with the following properties 

(AL-AH^h^^if = 7m.Q^'V, (108) 
iAL-AH^)P^,if = P^,Q(^V, 

for an arbitrary scalar (p. From (jl04p we obtain 

= K^{Al - AH')T,^ + - 3) ir±(Q(^))T - ^l,u^T 

4v,V.(ir,(Q(^)) + A) J_T. (109) 
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To get to the second line of ()109p we expressed T'^J in terms of T^^, and T using 
(jl(J6|) and the properties (|1U8|) of A^,. In the last term of the third line, the tachyonic 
pole Q(^) = could lead to a boundary divergence, impeding integration by parts 
when we contract with a conserved T^^,. 

Below we discuss the two cases of the (— ) and (+) solutions separately. In the 
(— ) case (this is the only case that has been so far discussed in the literature), 
there is the double pole present in the Green's function. Usually a double pole 
can be recast into a sum/ difference of two simple poles, giving rise to ghosts or 
tachyons. In the present case, this decomposition cannot be done so clearly. This is 
because nonlocal operators (such a square roots of the covariant d'Alembertian) are 
involved. The residue of the double pole diverges, which may or may not be a signal 
of a singularity of the perturbative approach at scales ~ H. On the other hand, the 
repulsive nature of the additional contributions (as compared to massive gravity) 
may be suggestive of a ghost-like state. However, given the results of Section 2, it 
is not entirely clear whether this ghost-like state, even if present, can be emitted in 
a final state, or it only appears as an intermediate state. In the latter case it would 
have been harmless. To clarify this, one needs a BRST invariant construction of 
the physical Hilbert space. Since the truncated theory considered here has only a 
limited interest anyway, we won't pursue this tedious project. One is certain, the 
conclusion reached in the literature that the amplitude has a ghost, is not entirely 
obvious beyond a reasonable doubt. 

For the (+) solution, the 1/ Q*-^^ divergence is not present since i^'+(Q*-^-') +A/H'^ 
also vanishes in the Q!^^'^ limit. Therefore, integrating on both sides of (jlU9p 
J dl^x^ T^^ ■ ■ ■ , for the (+) solution, (and integrating by parts the last term) 

The pole at Q*^^-* = has zero residue, and the residue of the pole at Q*^^^ = 
2H^ (^4 = —QH"^) together with the one at Ax, = give a massless graviton 
amplitude! Both, the tensorial structure in the numerator and, very importantly, 
the pole in the denominator, correspond to a propagating massless spin-2 state on 
dS background. What is this massless state? It corresponds to non-normalizable 
zero- mode discussed in Section 2. To understand this result better, it is instructive 
to look at the spectral representation of the amplitude . Let us rewrite 
as 

A= J T^^'K^ {Al - AH^) T^, - ^T'K+{Q)T - ^T'^K+{Q)T - It'-^T.(III) 
Which can be compared with the amplitude for massive gravity on de Sitter space 
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found in |52j: 

Am = j T'^"" {Al - 4:H^ - - 2H^)~^ Tf,^ - {Q + - 2H^y^ T - 

It'—(Q + M^-2H^)~^T ^T'-T. (112) 

It is sufficient at this point to obtain the spectral representation at y = which 
is completely determined by the function Kj^{z) for the ffist two terms in the r.h.s 
of (jllip . Below, z will refer to Q, and we will define the complex coordinate w = 
—H'^/A — z which in the flat space-time limit is positive in the physical domain, 
i.e., w — p2 > 0. Therefore, we study the spectral decomposition of K^{w) = 
K^(—w — -^) in the complex w plane. 

= m nlr. . . (113) 



—w — 



4 



(1 + |v^) 



where we have restored rric to be able to take the flat limit keeping rric finite. The 
fiat limit of this expression, K'^{w), gives the fiat propagator of the (+) branch 
described before: 

kl{w) = (114) 

This fiat propagator can be decomposed as a sum over KK contributions with spec- 
tral density p via a contour integral which goes around the cut just below and above 
the w positive real semi-axis and closes at infinity through a counterclockwise circle 
(the only contribution coming from the jump of the imaginary part of K across the 
cut, AlmK): 



+ ^ ^ Jo S-W^^' 271 Jo S-W 7T Jq s-ws^+m^^ ^ 

where (I114j) has to be interpreted as 

Kliw) = ^--= = ^ ^, (116) 

that implies that there is no pole for w = \w\e^'f' in the ffist Riemann sheet, i.e., 
< (j) < 2tt. The state corresponding to the pole of the propagator ()114p at 
w = e~'^'^m?^ is on the second Riemann sheet and does not contribute to physical 
amphtudes. Therefore, by a suitable prescription (which amounts to a choice of 
time direction) we are able to find an amplitude with no tachyonic instabilities. 

However, in trying to extend this argument to the curved background case, we 
face a problem. The curved generalization for pi3|l yields the following spectral 
representation: 

K+[W) = rH / Ti . 117 

Aw + H^ Ti Jo s-w{H{H + 2m^)+Asf + Amls 
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As in the flat case, the zero mode from the pole at w = e~'^^{H + 2mc)^/4 (corre- 
sponding to 2; = 2H'^ when rric = H) does not contribute in (|117|) while the KK's 
< wkk < 00 (— < zkk < 00) and the isolated mode w = e^H'^/A do. 
This prescription, also made the flat amplitude tachyon free, has turned the curved 
amplitude into the one of the (— ) branch with a ghostlike double pole in the third 
term of (ITTT|) . 

In order to describe the (+) branch amplitude, with only a massless pole an a 
branch cut of KK's, we must consider the other possible prescription, i.e., \/—w = 
VtT^^, the decomposition is then 



Now, the zero mode from the pole of pi3p at w = e*'^(if + 2mc)^/4 is on the flrst 
Riemann sheet and does contribute to the decomposition and the continuum enters 
with opposite sign. This relative sign implies that once the KK modes are treated as 
conventional positive-residue states, the residue of the zero mode has to be negative. 
This would imply that the zero-mode is a ghost. However, the zero-mode is not a 
normalizable state and because of this its interpretation is obscure. In particular, 
whether or not it can be emitted in a flnal state by any mater source localized on 
the brane is not clear. 

5 Outlook 

We studied the question of small perturbations on the self-accelerated solution of the 
DGP model. The issue is rather involved. We showed that the small perturbations 
on an empty self-accelerated background can be quantized without yielding ghosts 
or other instabilities. For conformal sources, such as radiation, small perturbations 
are instability free as well. More realistic, non-conformal sources, however, could 
trigger ghost-like instabilities in the linearized theory, and become unstable. We 
have suggested, in Section 4, possible loopholes in the latter conclusion, even within 
the linearized theory itself. More importantly, however, those non-conformal sources 
lead also to the breakdown of the linearized calculations, and, therefore, no conclu- 
sion on the (in) stability of the solution can be drawn without invoking non-linear 
dynamics. For the available non-linear spherically symmetric solution the dangerous 
linearized modes are spurious. 

To establish stability of generic non-conformal sources, one has to perform calcu- 
lations of 5D perturbations about a background produced by the source itself, since 
the latter is dominating over the cosmological background locally. One should see 
which of the obtained perturbations could be matched, at the scale r^, to solutions 
in the far-away region. Further interesting question is to study deviations from the 
self-accelerated background by introducing a small cosmological constant to explore 
the spectrum of the theory. We hope to report on some of these issues elsewhere. 
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Appendix A 

The goal of the present section is to discuss a toy model of a real scalar field min- 
imally coupled to a gravitational self-accelerated background. The idea is that the 
spectrum of small perturbations of this scalar bears some similarities to that of the 
gravitational perturbations. The action of the toy model is 

S = j d^xdy^l^-^g^^'dA^^dB^^^ + 2r, J d^'x^ (^-^r^d^ifid^ifi^ , (119) 

where A,B = 0, 1, 2, 3, 5, fj.,u = 0, 1, 2, 3,, (p{x) = y = 0), j^,i,{x) = g^u{x, y = 
0), and the metric of the spatially- flat self-accelerated background is 

ds'^ = (1 + H\y\)^ [-dt^ + a^{t)dx^) + dy"^ = A^{y)^^^{x)dx^dx'' + dy^ . (120) 

The equation of motion reads: 

dA (Vgg^^'dB^) + 2rJiy)d, {^^d,^) = . (121) 

Since ^/g = A'^a^ and ^^7 = ci^ we obtain for the bulk equation: 

-^d, iVir''d.<^) ^ □4$ = -M^M) . (122) 

Moreover, for the junctions condition we get 

[dy^]tl = -2r,n^^\o. (123) 



We look for the KK modes that are solutions to p22|) and ()123|) and can be expressed 
in the form: 

$(x,y)= ipm{x)fm{y)dm, □4V9m(x) = mVm(x) . (124) 



For the KK modes the bulk equation ()122|) and the junction condition ()123|) can be 
rewritten as follows: 



m'Uy) = -M^^ihl , (125) 
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and 



[dyfm{y)]t: = -2r,m^U0) • (126) 

As in the graviton case, and because of the boundary condition that depends on the 
mass of a KK mode ()126|) . the modes themselves are not orthogonal 

To solve the above equations it is useful to make the following change of variables 
from y to z, and introduce a new functions Um'- 

dz^^^, A-^l^u^ . (127) 

In terms of the new coordinate A{z) = exp{H\z\), and dy = A~^{z)dz- As a result, 
the bulk equation takes the form: 

^2^™^ /9^_^2\ ^^Q_ (128) 



dz^ 

The junctions condition p26|) rewritten in terms of Um reads: 

[d,uUz)]tl = {^H - 2r,m2) n^(0) . (129) 

From the above two equations it is immediately clear that there is a zero-mode 
solution m = 0, u = exp{3H\z\/2). However, this mode is diverging in the bulk. 

Other solutions of the above system of equations should be supplemented by 
the normalization condition, which in terms of the old and new variables reads as 
follows: 

dy^g''fl{y)= dzuUz) <oo. (130) 

Let us analyze the above system of equations. It is straightforward to check that 
there exists only one normalizable solution: 

3H 

^im. = e~"l^l , " = ^ "^c, (131) 

while the mass of this mode is 

ml = 3Hm^ - ml . (132) 

The above corresponds to the self-accelerated background if if = m^.. Then, we get 
a = H/2, ml = 2H^, and = {H/2Ay/\ 

What are the other solutions of ()128j) and ()129|1 ? The rest of the solutions 
constitute a continuum of massive states with > dH'^/A with the plane-wave 
normalizable wavef unctions: 



Um{z) = Asin{uj\z\) + Bcos{ujz) , u = \l m'^ —, (133) 
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where the normahzation coefficients A and B can be calculated 

A = -B ?!!^!^^ , (134) 

while the expression for B determines an important dependence of the wavefunctions 
on on the brane, in analogy with the flat case 



Appendix B 

In this appendix we consider linearized 4D massive gravity about a dS space with the 
special relation between the mass and dS curvature given in ()23|) . The Lagrangian 
density takes the form (see, e.g., jl2j, but the overall coefficient is different here) 

-\H'[K^ + \h^-\^l{hl^-h') . (135) 

For a generic value of the graviton mass the above Lagrangian has no gauge sym- 
metry. However, for a special value ()23p there is a symmetry present: 

5Ku{x) = (V^V, + ifV)p(a^) • (136) 

The above gauge transformation allows one to gauge-remove one out of five on-shell 
polarizations of the massive graviton However, this theory does not admit a 
coupling to non-conformal matter. There are many ways to see this. The most 
straightforward one is to write down the equations of motion, use the Bianchi iden- 
tities, and observe that the trace equation would necessarily imply = 0. This 
can also be established by noticing that introduction of the coupling h^^T^'^ in the 
Lagrangian p35j) violates the symmetry p3fj|l unless = 0. Then, one could in- 
troduce a coupling a la Stiickelberg (/i^j, — (V^V^ + if^7^i.)cr)T^'^, where a is a field 
that also transforms as 5a {x) = p{x) when the metric is transformed according to 
()136|] . Such a coupling would not violate the symmetry of the theory. However, the 
equation of motion of the a field would require that = 0. As we have seen in the 
Section 2, in the DGP model this difficulty is avoided due to the extrinsic curvature 
terms. 

To this end we would like to calculate the propagator of the theory given by (jl35|) 
without coupling it to any matter. The gauge independent part of the propagator 
would tell us how many degrees of freedom are propagating in the theory. Without 
introducing the source, the propagator is defined as an inverse of the quadratic 
operator in p35p . However, because of the symmetry ()136|) . this operator is not 
invertible. Hence, we should proceed by introducing the gauge fixing term which we 
choose to be 

C„f = --h^ . (137) 

a 
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After this term is included, we introduce a "spectator" T^jy, for which T 7^ 0. For- 
mally, this is possible because of the gauge fixing term. However we should keep 
in mind that the original theory wouldn't admit any non-conformal matter, and, 
in that respect, the introduction of the "spectator" non-conformal source after the 
gauge fixing is just a technical trick that enables one to extract the non-derivative 
terms of the propagator without too much of calculations. We follow this philosophy 
below. 

In this gauge the equations of motion read: 

( \ \ 1 

-H'^ h^^ + -l^yh - — i {h^^ - 7^^/i) + --i^yh = -T^^ . (138) 
\ 2 / Z a 

The Bianchi identities enforce the following constraint: 

777^ 1 

^(V^ V - V,/i) = —V,h . (139) 
2 a 

The trace of equation ()138p takes the form 

-H^\h + -h= -T . (140) 

The above system of equations can be solved. In what follows we concentrate on 
the case when ()23|) is satisfied, i.e., ml = 2H^ . First we solve for the trace 

h = -. (141) 

As we see, in a theory without the gauge fixing term (i.e, when a —>■ 00) the field 
diverges. 

Next using the identity we solve for h^^: 

Fortunately, the double pole terms on the r.h.s. are gauge dependent On the other 
hand, the gauge-independent part of the above expression gives the non-derivative 
part of the inverse of the gauge kinetic term of ()135p . The latter takes the form 

2 ^4 - ^ ' 

The tensorial structure in the numerator of Eq. ()143|) confirms that the theory p35|) 
at the linearized level does not propagate the fifth, helicity-0, degree of freedom 



^^Notice that the gauge-dependent terms in this case contain not only derivatives but also terms 
proportional to 7^,^. This is a direct consequence of the form of the gauge transformations H136|l . 
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Appendix C 

In this appendix we give the expressions used in Section 3 for the first and second 
order perturbations of the Ricci tensor around a de Sitter background metric 7^1,. 
The metric is expanded as 

9,,u = 7mv + 5g^^ , (144) 

with inverse 

g,u^Y'_Sg^'' + 5g^^5g'^\ (145) 
The first order Ricci tensor components are given by 

RilKSg) = -\ - V«V^5(7° - V^VJg^^ + V.VJgZ) , (146) 

which for a transverse traceless perturbation, 5g^y = hJ^J, reduces to 

Ri^ih^^) = -lahl^ + ^HXI ■ (147) 

Note that the Vq, has to be commuted through V^^i, in the right hand side of p46|l 
to get ()147p . Here, we recall that for the de Sitter background metric 7^1, we have 

The second order terms are given by 

R?^ (Sg) = +\ i^'Sg: - 2V"SgPJ {\/,5g,p + WJg.p - V,Sg,,) 

+lsg'"'{V^VJgpu + V^VJgpa-V^Vp5gua) • (148) 

Which for a transverse traceless perturbation reads (we drop the TT superscript on 
the r.h.s.) 

-\h'"' {yaV^.hp, + V„V,V - VaVp/i^, - V^Nahp, - V^V./ip^ + V pKa) (-149) 
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